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METHOD AND SYSTEM FOR WAVELET-BASED 
REPRESENTATION AND USE OF CARTOGRAPHIC DATA 



TECHNICAL FIELD OF THE INVENTION 
The present invention relates generally to cartography, and in particular, to the use of 
wavelet transforms to represent cartographic data. 

BACKGROUND OF THE INVENTION 
Cartographic databases are used in computer-based systems to represent geographic 
features on displayable maps. These geographic features may include roads, intersections, and 
points-of-interest, as well as other features, such as lakes, railroad tracks, buildings, airports, 
stadiums, parks, mountain ranges, bridges, and so on. To permit these features to be displayed, 
cartographic databases usually include geographic location information in a format that allows 
computers to display maps illustrating the size, shape, position, and character of these various 
geographic features on a display device. 

In order to represent these continuous features on a computer, it is necessary utilize some 
form of discretization. Typically, the features are represented in terms of discrete data points 
consisting of coordinate pairs (longitude, latitude) or coordinate triples (longitude, latitude, 
altitude). A point feature such as a point of interest or intersection is then represented as a single 
data point, while other features such as roads, lake boundaries, and so on are represented as a 
sequence of data points. The endpoints of the geographic features are commonly referred to as 
nodes while data points between the nodes are referred to as shape points. Accordingly, nodes 
and shape points can be used to approximate the true shape of the geographic feature. The shape 
points and nodes are stored in what is known as a shape point database, which may be part of a 
larger cartographic database. With this representation scheme, the spacing between shape points 
is selected such that the error between the data point representation and the actual geographic 
feature lies within some acceptable bound. In order to achieve a high level of accuracy and 
detail in the data point representation, it is often necessary to store a large number of data points. 
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A function available with some computerized map display systems is the ability to zoom. 
Zooming is the process of magnifying or reducing the scale of a map displayed on the monitor to 
change the level of detail. Reducing the display scale in a displayed map is referred to as 
zooming out, while increasing the display scale is referred to as zooming in. 

To permit zooming, some cartographic databases include separate collections of location 
data, each with a different level of detail, which can be used by the map display function. Using 
these different layers of cartographic data, the map display function can provide or support 
zooming. However, providing multiple layers of cartographic data requires additional database 
storage space. 

In one type of display system, the different resolution layers of cartographic data are 
produced by using line generalization algorithms (also referred to as line simplification 
algorithms). Many different line generalization algorithms are currently being used in 
cartography, one of which is the Douglas-Peucker algorithm. Generally, line generalization 
refers to a process of reducing the data points representing a linear geographic feature in a map, 
while maintaining the overall shape and/or trend of the feature. By applying a line generalization 
algorithm to a cartographic database having a large number of data points representing a high 
level of detail, lower resolution databases can be generated. Although known line generalization 
algorithms are adequate for producing cartographic data collections with different levels of 
detail, there is an opportunity to improve upon the speed and efficiency at which the algorithms 
operate. 

Accordingly, there is a need for a method and system for representing cartographic data 
that improves upon known generalization and zooming techniques and provides an ability to 
reduce the storage requirements of the data without sacrificing feature resolution. 

BRIEF DESCRIPTION OF THE DRAWINGS 
FIG. 1 is a Cartesian graph showing a Haar basis function and scaling function. 
FIG. 2 is series of Cartesian graphs illustrating an exemplary square integrable function 

f(x) and its approximation at successively higher resolutions using wavelet basis functions of 

successively higher resolutions. 



FIG. 3 shows a pair of graphs illustrating an exemplary wavelet basis function and the 
corresponding wavelet scaling function. 

FIG. 4 is a block diagram illustrating an exemplary system for converting shape point 
data of a cartographic database into a database of wavelet and scaling coefficients in accordance 
with an embodiment of the present invention. 

FIG. 5 is a flow chart illustrating a method of operating the processor shown in FIG. 4. 

FIG. 6 is a block diagram illustrating a computer-based system capable of displaying 
cartographic data represented using wavelets in accordance with a further embodiment of the 
present invention. 

FIG. 7 is a flow chart illustrating a method of operating the computer-based system 
shown in FIG. 6. 

FIG. 8 is a block diagram illustrating a system for computing multi-resolution error 
metrics in accordance with another embodiment of the present invention. 

FIG. 9. is a flow chart illustrating a method of operating the processor shown in FIG. 8. 

DETAILED DESCRIPTION OF THE PRESENTLY 
PREFERRED EMBODIMENT(S) 

To address the needs for reducing cartographic database sizes and, improving upon 
known line generalization and zooming techniques, the present invention provides a method and 
system for generating and using wavelet representations of cartographic data. 

The wavelet transform maps a square integrable function to the scale-space domain. 
Specifically, the wavelet transform of a signal provides information about the resolution content 
of a function as well as information about the location of each resolution component. This 
property of wavelets is called time-frequency localization. Wavelet transforms are widely used 
for the analysis of non-stationary signals. Their unique properties make them a valuable tool in a 
wide variety of areas such as numerical analysis, edge detection in images, image compression, 
and analysis of fractal signals, to name just a few. 



A wavelet representation of cartographic data can provide an accurate representation of 
the data, and a consistent and intuitive method for generalization for the purpose of map display 
at different zoom levels. Furthermore, the wavelet transform can be easily extended to two or 
5 more dimensions. 

Display operations like generalization and zooming can be accomplished much more 
intuitively and easily using the wavelet representation, because generalization will simply 
involve dropping wavelet coefficients at certain scales. With a wavelet representation, the 
wavelet coefficients for each geographic feature can be stored for the highest level of resolution. 
10 To produce displays of lower resolution, higher resolution wavelet coefficients are not used. 

•ji Wavelets and Wavelet Transforms 

J; Wavelets are families of functions of the following form 

IU / \ i 1-1/2 (x-b\ /1X 

y/ ab (x) = \a\ y/ (l) 

s V a ) 

ft 5 where a and b are respectively called the dilation (scaling) and translation parameters. As a is 
increased, the function y/ ab (x)becomes more spread out or dilated, and as b is increased or 

O decreased y/ ab (x)is shifted or translated by b I a . Thus, wavelets are constructed by dilating and 

translating a single prototype function y/{x) called a mother wavelet. The dilation and translation 

parameter pairs (a, b) can range over a continuous or discrete subset of x 5R . 

20 If the parameters a and b are restricted such that a = 2~ m and b = mi , where m and i 

are integers, then the resulting wavelets are called dyadic wavelets and are parameterized by the 
integers m and i . Dyadic wavelets are expressed as: 
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^W = 2 m/ V(2 m x-0 (2) 

In this case, negative values of m correspond to dilating y/{x) and positive values of m 
correspond to contracting y/{x) . 

Although the present invention is not so limited, the embodiments described herein can 
utilize dyadic wavelets. It is possible to select a dyadic mother wavelet y/(x) such that the set of 
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functions {i// mi }v(w, i) e Z x Z constitutes an orthonoraial basis for L 2 (9t) , where 



L 2 (9t) represents the space of square-integrable functions defined on the real line. The simplest 
choice of y/(x) which yields an orthonormal basis of Z 2 ($R) is the Haar mother wavelet, shown 
in the FIG. 1. 

Other orthogonal bases include the Daubechies family of bases. These functions are 
indexed by an integer N ; the regularity (which loosely means smoothness) of these bases 
increases with N. 

Turning now to wavelet transforms, which rely on wavelets such as those discussed 
above, there are many different types of wavelet transforms. The wavelet representation of 
geographic features yields accurate results, and can also be made more accurate by using 
smoother basis functions, and using more coefficients in certain regions of each geographic 
feature. 

The different types of wavelet transforms are distinguished from one another by the types 
of wavelets used in their definitions. In the methods and systems disclosed herein, orthonormal 
wavelet transforms are preferably used; however, nonorthonormal wavelets could also be used 
for the purpose of representing cartographic data. The semi-discrete orthonormal wavelet 
transform of a square integrable computable function f(x) is defined, using dyadic wavelets, by 
the analysis and synthesis equations given below: 

Analysis Equation: 




(3) 



-oo 



Synthesis Equation: 



00 00 



/(*)= I 2n™,i)y mi (x) 



(4) 



The values F(m,i) are called the wavelet coefficients. The functions 
Vmi i x ) = 2 m/2 ^(2 m x - i) are the wavelet basis functions. The set of basis functions \t// mi (x)} is 
thus constructed by scaling and translating the mother wavelet function y/{x) * The integer m is 
the dilation or scaling parameter and the integer i is a translation parameter. Negative values of 
5 m correspond to dilating y/(x) and positive values of m correspond to contracting ^(x) . 

As can be seen from the analysis equation above, the wavelet coefficients F(mj) are 
obtained via an inner product computation of the square integrable function f{x) and the 
corresponding wavelet basis function, i.e., (f(x\ y/ mi (x)) . The set of wavelet coefficients 
Q {F(mj)} can be indexed by both scale or resolution m and location or translation i . 
if The square integrable function f(x) can be reconstructed using the synthesis equation, 

I?. which is sometimes called a wavelet series expansion. The wavelet representation of a square 
vy integrable function is thus a weighted sum of wavelet basis functions, the weights being the 
m wavelet coefficients. 

j|\ The wavelet transform is considered to be semi-discrete because it takes the continuous 

j| function f(x) and transforms it into a set of coefficients indexed by integral values. This 

transform therefore accomplishes the goal of representing the continuous function f(x) in terms 
O of a set of discrete values {F(m, /)} , which can be stored in memory on a computer. 

The limits on the translation index i in the inner summation of Equation (4) are 

determined by the support of the square integrable function f{x) . The limits on i are finite if 
20 f{x) has compact support. The limits on the dilation index m in the outer summation of 

Equation (4) are infinite. In practice, however, the limits of m must be truncated both above and 

below to predetermined values to yield a reasonable approximation of f(x) . This is done by 

collecting all terms in the expansion for m e(- oo ? M') into a single term A M ,(x), which 
represents the approximation of the function at the highest resolution level, layer 0. The finer 
25 details of the function are expressed in terms of a wavelet series expansion as indicated in 
Equation (5). 
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M oo 

f(x)*A M .(x)+ y S^ko^-W 

m-M' /=-oo 



(5) 



It is often convenient to include only nonnegative scales in the approximation, i.e., 
truncate all expansion terms below m = 0 and above some arbitrary scale M that yields an 
5 approximation to f(x) that is acceptable for the given application. Furthermore, A M ,(jc)is 

expressed in terms of dilations and translations of the function y/{x) , which is called the scaling 
function associated with the wavelet function <j>{x) . This results in the following equation: 

co M oo 

a /to* Z^(OA«to+ £ I^("M> m ,to (6) 

""„[ i'=-oo m=M' i=-co 

fi) where the quantities the F'(z') , are known as the scaling coefficients, and are given by 



*"(/)= J/(x^(x>fc (7) 



3?f One advantage of using the wavelet transform for representing cartographic data is that 

U the transform organizes the information in the data into components at various resolutions or 

scales. This is particularly useful for line generalization that supports map display features, such 
as zooming. This can be explained as follows. For negative values of m , the corresponding 
wavelet basis functions are dilated, i.e., they have a larger width of support than the mother 
20 wavelet. As a result, the wavelet coefficients indicated by negative values of m will correspond 
to coarse scale or low-resolution information contained in the cartographic data. This 
corresponds roughly to zooming out and observing the square integrable function f(x) from a 
distance. The more negative the value of m , the more the displayed feature is effectively 
zoomed out. 
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On the other hand, for positive values of m , the resulting wavelet basis functions are 
contracted with respect to the mother wavelet, i.e., their width of support is smaller than that of 
the mother wavelet. As a result, the wavelet coefficients indicated by positive values of m 
correspond to fine scale or high resolution information contained in the cartographic data 
represented by the square integrable function f(x) . This corresponds roughly to zooming in and 
observing the square integrable function f(x) from close up. 

This zooming aspect of wavelet transforms is illustrated in FIG. 2 using the Haar basis 
function. The column of graphs on the right-hand side of FIG. 2 shows an exemplary square 
integrable function f(x) and its approximation at successively higher resolutions using wavelet 
basis functions of successively higher resolutions. The corresponding wavelet basis functions 
for the successive levels of resolution are illustrated in the left-hand column of graphs. 

From the sequence of graphs, it is apparent that as each successive scale of wavelet 
coefficients is added to the approximation, the approximation becomes more accurate. The Haar 
basis function is used here only for illustrative purposes. There exist other wavelet bases that 
can yield accurate representation of continuous functions. 

FIG. 3 shows a pair of graphs illustrating an exemplary wavelet basis function and an 
exemplary wavelet scaling function. The wavelet basis function shown in FIG. 3 is a 
Daubechies sixth order basis function. This is an orthonormal basis function that has compact 
support and has a high enough order of continuity for accurately representing cartographic data 
for some applications. The wavelet scaling function shown in FIG. 3 is associated with the 
Daubechies sixth order basis function. The Haar basis function is a degenerate case of 
Daubechies basis function. Daubechies bases do not have closed form expressions, but since all 
wavelet basis functions are derived from a single prototype function, it is sufficient to compute 
the values of the mother wavelet once and then store them in a look-up table. There also exist 
other (non-compactly supported) orthonormal wavelet bases, such as the Battle-Lemarie basis 
function, which can also be used in other embodiments of the invention. 
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Other wavelet bases include compactly supported non-orthogonal wavelet bases, such as 
the B-wavelets which are constructed from splines. The Haar basis function is a degree zero 
B-spline. These may be of special interest in cartographic applications; however, since they are 
not orthogonal, the wavelet coefficients will not be quite as easy to compute, as compared to the 
orthonormal case. However, this is generally not an issue since the wavelet coefficients for a 
cartographic database need only be computed once, at the time the wavelet database is generated. 

Parameterization of Cartographic Data 
Each geographic feature can, in general, be represented as a continuous function in space 

and can be represented in implicit form as f(x 9 y 9 z) = 0 , where x , y and z represent longitude, 

latitude and altitude, respectively. In order to apply the wavelet transform, it is convenient to 

express each geographic feature in parametric form as a group p(*) of polynomial functions of 

an independent parameter t as follows: 

p(0=M0 y(t) *(')] (8) 

The parameterization scheme can be chosen in such a way that each of the functions x(t) , y(t) 
and z{t) are each single valued. The wavelet transform can then be applied to each of the 
functions x(t) 9 y(t) and z(t) . This yields P(mj)=[x(mj) Y(mj) z(m,i)] where X(m 9 i) 9 
Y(m 9 i) and Z{m 9 i) are the wavelet coefficients of x(t) , y(t) and z(t) respectively, and 

P^0=[jif'(0 7 (0 z (0] where ^ (0> ^'(0 and z (0 ^ e the scaling coefficients of x(^) , 
y(i) and z(t) , respectively. 

System Architecture 

FIG. 4 illustrates an exemplary system 20 for converting a shape point database into a 
database of wavelet and scaling coefficients in accordance with an embodiment of the present 
invention. The shape point database itself may be part of a larger cartographic database, which 
may include other kinds of databases used for functions such as route calculation, maneuver 
generation, and other navigational functions. The system 20 includes a processor 22 in 
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communication with a shape point database (DB) 24 and a wavelet database 26. A user interface 
21 permits database developers to configure conversion parameters used by the processor 22 
such as the desired highest level of resolution in the wavelet representation. 

The shape point database 24 stores a plurality of shape points and nodes representing one 
or more geographic features that can be displayed by a computer-based system, such as a 
navigation system, as part of a digital map. Shape points and nodes represent geographic 
locations within a map. In the embodiments disclosed herein, the shape points and nodes can 
consist of coordinate triples (longitude, latitude, altitude) or coordinate pairs (longitude, latitude), 
which represents specific locations along geographic features. Nodes are special cases of shape 
points, representing, for example, intersections between road segments and end points of 
geographic features. 

Data stored in the shape point database 24 can be organized on a per geographic feature 
basis, i.e., each geographic feature is represented as sequence of nodes and shape points. 
Although the database 24 can store a number of geographic feature entities 30, three entities 
SEG 0 , SEGi, SEG 2 are shown for illustration purposes. A pair of nodes and zero or more shape 
point entries 31 can be associated with each entity. For entities representing straight or 
substantially straight geographic features, the nodes 33 are sufficient to accurately represent the 
feature; and therefore, for straight geographic features, the nodes alone can be stored without 
shape points. 

In the example shown, a first geographic feature SEG 0 is represented by shape points 
SPo(O) - SP 0 (i), where i is an integer, and nodes EN 0 (0) - EN 0 (1); a second feature SEGi is 
represented by shape points SPi(O) - SPi(j), where j is an integer, and nodes ENi(O) - ENi(l); 
and a third feature SEG 2 is represented by shape points SP 2 (1) - SP 2 (k), where k is an integer, 
and nodes EN 2 (0) - EN 2 (1). 

The processor 22 can be any computing device for performing the steps shown in the 
flow chart 50 of FIG. 5, including a conventional personal computer (PC) executing a software 
program. In such an embodiment, the software program can be stored on a computer-usable 
medium for configuring the processor 22 to perform the method of FIG. 5 as disclosed herein. 
The computer-usable medium can be any article of manufacture for storing computer programs, 
such as a CD ROM, hard drive, DVD, diskette, solid-state memory such as a RAM, ROM, or 
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EEPROM, or the like. Alternatively, the processor 22 can be a custom hardware device, such as 
one or more application specific integrated circuits (ASIC) designed to perform according to the 
method shown in FIG. 5. 

As will be discussed in greater detail below, the processor 22 processes data stored in the 
shape point database 24 to convert it into a plurality of wavelet coefficients that are stored in the 
wavelet database 26. The wavelet representation can provide a more efficient representation of 
map data that facilitates line generalization and map display features, such as zooming. 

Generally, for each geographic feature entity in the shape point database, the processor 
22 generates a corresponding set of wavelet and scaling coefficients that are stored in the wavelet 
database 26. The geographic features can be reconstructed using the wavelet synthesis equation 
to draw the features on a map display. Further, the wavelet and scaling coefficients for each 
entity can be indexed by layer or resolution level. 

In the example shown, the first geographic feature SEGo is represented by wavelet 

coefficients P 0 (0,0), . . . , P 0 (M 0 ,N q ) 9 and scaling coefficients P 0 ' (o), . . . , P 0 (iV 0 ), where M 0 and 
N 0 are integers; the second geographic feature SEGi is represented by wavelet coefficients 

P 1 (0,0),...,P 1 (M 1 ,A^ 1 ) and scaling coefficients ^ (o),...,^^^), where M x and N } are integers; 
and the third geographic feature SEG 2 is represented by wavelet and scaling coefficients 
P 2 (0), . . . , P 2 (iV 2 ) , where M 2 and N 2 are integers. In the above, 

P k (m 9 i)=[X k (m 9 i) Y k {mj) Z k (mj)] and P k (i)= [x^(i) 7/(z) Z 4 '(i)]. (9) 



The shape point and wavelet databases 24 and 26, as well as all other databases disclosed 
herein, can be stored using any suitable computer-usable storage media, such as one or more CD 
ROMs, hard drives, DVDs, diskettes, solid-state memories such as RAMs, ROMs, or 
EEPROMs, or the like. 

FIG. 5 shows a flow chart 50 illustrating a method of operating the processor 22 shown in 
FIG. 4. In step 52, the processor 22 receives shape points and/or nodes representing one or more 
geographic features from the shape point database 24. The processor can include a memory (not 
shown) for buffering the shape points. 
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Next in step 56, the wavelet coefficients are computed for each parameterized geographic 
feature p(t) = [x(t) y(t) z(t)]. This can be done in one of two possible ways. If p(*) is known 
at all points, or if the data points of p(r)are spaced sufficiently closely, then the wavelet and 
scaling coefficients are computed for the given data points by using the following discretized 
versions of Equations (3) and (7): 



F(m,i)* ^f{k^)y/ mi {kl^x (10) 



Jt = -00 



tsfikAx^JkAx)^ (11) 



£=-00 



If the data points of p(t) are not spaced sufficiently closely, then the wavelet coefficients 
can be computed by using a least squares fitting technique. 

Many different types of wavelet transforms are available and the present invention is not 
limited to any particular type. However, a preferable form of the wavelet transform employed by 
the method shown in FIG. 5 was described above in connection with Equations (3)-(4). 

An advantage of the wavelet representation of a function is that the wavelet and scaling 
coefficient set often contains a large number of small or zero values, which can simply be 
dropped. This allows cartographic data being represented by wavelets to be compressed, 
reducing the amount of database space necessary to store the wavelet coefficients, without 
significantly degrading the accuracy of the representation. 

In step 58, the wavelet coefficients and scaling coefficients for the geographic features 
are stored in the wavelet database 26. 

FIG. 6 illustrates a computer-based system 100 capable of displaying cartographic data 
represented using wavelets in accordance with another embodiment of the present invention. 
The system 100 includes a computer 102, a wavelet database (DB) 104, a user interface 108, and 
an output display device 106. 
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The computer 102 can be any computing device, including a conventional PC or other 
microprocessor-based device having a processor 1 12 and a memory 1 10. The processor 1 12 can 
execute a computer software program stored in the memory 110. The software program can 
direct the processor 1 12 to perform the steps shown in the flow chart diagram 120 of FIG. 7. 

The user interface 108 can be a software-generated graphical interface included in a 
commercially-available PC operating system, configured to direct the processor 1 12 to receive 
and process user selections, such as the desired level of resolution in the displayed map. 

The wavelet database 104 includes the wavelet coefficients and scaling coefficients and 
basis function data. 

The computer-based system 100 of FIG. 6 can be part of a navigation system included 
with on-road vehicles, such as automobiles, trucks, buses, and so on, as well as advanced driver 
assistance systems (ADAS) provided in on-road vehicles. 

The flow chart diagram 120 in FIG. 7 illustrates the operation of the computer-based 
system 100 shown in FIG. 6. In step 122, the computer 102 receives a user selected resolution 
level from the user interface 108. The resolution level determines the scale at which the 
cartographic data is displayed, and may be embodied as a zooming feature of the system 100. 
The zooming feature is described in further detail below. 

Next, in step 124, the processor 1 12 retrieves wavelet and scaling coefficients 
corresponding to the desired geographic features from the wavelet database 104. The wavelet 
coefficients in the database 104 can be indexed by scale, facilitating retrieval of the coefficients 
corresponding to the user selected resolution level, if desired. In addition to retrieving the 
wavelet coefficients, the processor can also retrieve from the database 54 the precomputed 
wavelet basis and scaling function values. The parameterized functions 

p(f ) = [x(t) y(t) z(t)] , which represent geographic features, such as a road segments, rivers, 
or the like, for cartographical applications, can be reconstructed from the wavelet coefficients 
and a look-up table storing the precomputed wavelet basis and scaling function values. 
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In step 126, the processor 112 reconstructs the functions that represent geographic 
features. Through a conventional display interface card (not shown), the processor 112 then 
passes the reconstructed functions to the output display device 106, which then displays the 
function representing the geographic feature to the user (step 128). The wavelet zooming feature 
5 of the system 100 is explained as follows. Starting with the highest resolution wavelet series 

approximation of the function p(t) = [x(t) y(t) z{t\ zooming out, or reducing resolution, can 
be accomplished by successively removing the finest scale wavelet coefficients from the wavelet 
expansion. As an example, suppose that the cartographic display data is organized into five 
layers, labeled layer 0 through layer 4, with layer 0 representing the lowest level of resolution 

lp and layer 4 representing the highest level of resolution. Let p.(/) denote the function 

representing the approximation of the function p(^) resolution i . Then, p 4 (t) is obtained by 
y including wavelet coefficients for scales 0 < m < 4, i.e., f 4 (x) is a summation of the form of 
Vi 1 Equation (5), with M = 4 , as shown below in Equation (12). 

Yf'. Layer 3 data can be represented by dropping all wavelet coefficients corresponding to the 

%5 finest scale used in the initial approximation p 4 (t) , resulting in the wavelet series shown in 
Equation (13); layer 2 data can be represented by dropping all wavelet coefficients 

In corresponding to the finest scale in the layer 3 representation, resulting in the wavelet series of 
Equation (14); and so on. As scales are dropped from the wavelet approximation, the resolution 
of the approximation becomes progressively lower, i.e., the wavelet representation is effectively 

20 zooming out. Equations (12)-(16) show the wavelet synthesis equations for each layer. 



Layer 0: 


PoW* 


|;po> 0( .(x)+x|;p(m,o^(0 

t =-co m=0 i=-co 


(12) 


Layer 1: 


Pi(0* 


/=— co m=0 i=-<x> 


(13) 


Layer 2: 






(14) 


Layer 3: 


PaW- 


; Xp'(/>o,W+XSpko^ ; (0 

i=-oo m=0 *=-co 


(15) 


Layer 4: 


P4W* 




(16) 



i=-co 
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Multi-resolution Error Metrics for Cartographic Databases 
In addition to generating displayable maps, the wavelet coefficients of the cartographic 
data can be used to define a multi-resolution error metric for measuring database accuracy. 
5 FIG. 8 illustrates a system 200 for generating a multi-resolution error metric in accordance with 
yet another embodiment of the present invention. The database which is being assessed for 
accuracy is commonly referred to as a sampled database 206. The system also includes another 
database called the ground truth database (DB) 204 and a processor 202. Both of these databases 
are shape point databases, i.e., both the sampled database 206 and the ground truth database 204 
*H) contain plural shape points and nodes representing the same geographic features, such as roads, 
^ and so forth; however, the shape points are more densely represented in the ground truth 
O database 204 to give a highly accurate representation of the geographic features. The ground 
truth database is constructed by using more accurate sources (which may include preexisting 
^ maps or imagery, as well as sensor measurements) than those sources used for constructing the 
j 5 sampled database. 

12 The processor 202 can be any computing device, including a conventional PC or other 

G microprocessor-based device. The processor 202 can execute a computer software program for 
□ performing the steps shown in the flow chart diagram 230 of FIG. 9. 
"** The flow chart diagram 230 illustrates a method of operating the processor 202 to 

20 generate one or more multi-resolution database error metrics. In step 232, the wavelet and 
scaling coefficients of the geographic features represented in the ground truth database are 
computed and stored in the ground truth wavelet database 204. This can be done by generating 
the coefficients from the shape points of the ground truth database (GD) on a per segment basis 
as described above in connection with FIGS. 1-5. 
25 The wavelet and scaling coefficients of the geographic features represented in the 

sampled database are then calculated and stored in the sampled wavelet database using the 
method as described above in connection with FIGS. 1-5 (step 236). 

The total error of the sampled database 206 with respect to the ground truth database 204 

is designated e(0 = [e x (t) e y (t) e z (t% where e x (t), e y (t) and e z (t) are the components of 
30 error corresponding to longitude, latitude and altitude, respectively. 
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In step 238, the total error e(t) is computed as follows; 

CO M CO 

•W=ZP^(0-P®(0^,W+ Z l^s D (rn,i)-P GD (m,i)\w mi it) 

l=-eo m =M '['=-« (17) 

= ®«"W+y" „.e.(t) 

where 

cO 

«„(') = Ef p ffl)( w '0-P OT ('«»0K ( W fOT M' <m<M (18) 

/=-co 

and P GZ) and are the wavelet coefficients of the ground truth database 204 and the sampled 
database 206, respectively and P GD and are the scaling coefficients of the ground truth 
database 204 and the sampled database 206, respectively. The set of values 

{ e m (0 = [ e m (0 e m (0 e m (*)]} represent the layer error metric values. 

In step 240, a layer error metric value can be generated for a selected resolution layer of 
the sampled database 206. The total error determined by Equation (14) is thus decomposed into 
error components at different scales, represented by the integer m. Thus two types of error can 
be obtained for each geographic feature represented in the sampled database: the layer error 
e m (t) at a particular scale m , as well as the total error e(t) , which is the sum of the errors at all 
scales. Note that B m (t) and e(t) are vectors, with longitude, latitude and altitude components. 
But we can compute e c m om bined (t) = (e x m (tf + e y m (tf + e z m (tf ) and 

e combmed ^ _ + e y ^y + ^ ^y j fo ^ a com ^ no ^ measure c f error i n a n three directions. 

In step 242, the values of the error metric are stored by the processor 202 for later 
analysis and use by a computer system. 



16 



Conclusion 

Wavelets are useful for yielding a compact representation of cartographic data with 
improved accuracy, and for line generalization, and zooming. Furthermore, different wavelet 
bases can be applied to a cartographic database, once the initial wavelet formulation has been 
done. Moreover, the wavelet transform can be extended to two or more dimensions according to 
teachings disclosed herein. For instance, a two-dimensional wavelet representation of 
geographical information can be used to represent surfaces. 

While specific embodiments of the present invention have been shown and described, it 
will be apparent to those skilled in the art that the disclosed invention may be modified in 
numerous ways and may assume many embodiments other than those specifically set out and 
described above. Accordingly, the scope of the invention is indicated in the appended claims, 
and all changes that come within the meaning and range of equivalents are intended to be 
embraced therein. 
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WHAT IS CLAIMED IS: 



1 1 . A method for representing cartographic data in a computer-based system, 

2 comprising: 

3 computing a plurality of wavelet and scaling coefficients corresponding to at 

4 least one function representing a geographic feature in a cartographic database; and 

5 storing the wavelet and scaling coefficients in a computer-usable database, the 

6 coefficients being usable for representing the cartographic data in the computer-based system. 

7 

3 2. The method of claim 1, wherein the geographic feature is originally represented 

*2 by a plurality of data points. 

4 

[J 3. The method of claim 2, wherein the data points are selected from the group 

S consisting of coordinate pairs and a coordinate triples. 

3 

d 4. The method of claim 1 , wherein the geographic feature is the boundary of a 

;| feature selected from the group consisting of a road, waterway, building, park, lake, railroad 

3 track, and airport. 
4 

1 5. The method of claim 2, wherein the step of computing the wavelet coefficients 

2 and scaling coefficients includes applying a wavelet transform to a function defined by the data 

3 points representing the geographic feature. 
4 

1 6. The method of claim 1 , wherein the step of computing the wavelet coefficients 

2 and scaling coefficients includes: 

3 computing the wavelet coefficients by performing a least-squares fit. 

4 

1 7. The method of claim 1 , wherein the wavelet and scaling coefficients are 

2 computed using a semi-discrete orthonormal wavelet transform. 
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1 8. A method of displaying on a computer output device a function representing a 

2 geographic feature, comprising: 

3 retrieving from a computer-usable database a plurality of wavelet and scaling 

4 coefficients associated with the geographic feature, the coefficients being derived from a 

5 plurality of data points specifying geographic locations according to a predetermined reference 

6 system; 

7 computing the function using the retrieved wavelet and scaling coefficients; and 

8 displaying the function on the computer output device. 

9 

1 9. The method of claim 8, wherein the data points are selected from the group 

<jf consisting of coordinate pairs and a coordinate triples. 

Eft 10. The method of claim 8, wherein the geographic feature is selected from the group 

f[? consisting of a road, waterway, building, park, lake, railroad track, and airport. 

H 1 1 . A system for displaying on a computer output device a representation of a 

i™2 geographic feature, comprising: 

^| a database storing a plurality of wavelet and scaling coefficients associated with 

13 the geographic feature, the wavelet and scaling coefficients being derived from a plurality of 

5 data points specifying geographic locations according to a predetermined reference system; 

6 a processor configured to calculate a function using the wavelet and scaling 

7 coefficients, the function representing the geographic feature; and 

8 a display device for displaying the function. 
9 

1 12. The system of claim 1 1, wherein the data points are selected from a group 

2 consisting of coordinate pairs and coordinate triples. 

3 
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1 1 3. A method of generating a computer-usable database that represents cartographic 

2 data using a plurality of wavelet and scaling coefficients, comprising: 

3 providing a predetermined database that represents the cartographic data using a 

4 plurality of data points specifying geographic locations; 

5 computing a plurality of wavelet and scaling coefficients from the data points; and 

6 storing the wavelet and scaling coefficients in the computer-usable database. 

7 

1 14. The method of claim 13, wherein the data points are selected from the group 

2 consisting of coordinate pairs and coordinate triples. 

J. 

J3 1 5 . The method of claim 1 3 , wherein the geographic feature is the boundary of a 

f4 feature selected from the group consisting of a road, waterway, building, park, lake, railroad 

0| track and airport. 

if! 

1 16. A system for generating a computer-usable database that represents cartographic 

^2 data using a plurality of wavelet and scaling coefficients, comprising: 

r$ a first computer-usable database storing the cartographic data represented using a 

JS| plurality of data points specifying geographic locations; 

13 a processor configured to compute a plurality of wavelet and scaling coefficients 

6 from the data points; and 

7 a second computer-usable database, operatively coupled to the processor, for 

8 storing the wavelet and scaling coefficients. 
9 

1 17. The system of claim 1 6, wherein the data points are selected from the group 

2 consisting of coordinate triples and coordinate pairs. 

3 

1 18. The system of claim 1 6, wherein the wavelet coefficients and scaling coefficients 

2 are computed by applying a wavelet transform to a function defined by the data points 

3 representing a geographic feature. 
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1 19. The system of claim 16, wherein the wavelet coefficients are computed by 

2 performing a least-squares fit. 

3 

1 20. A method for generating a database error metric in a computer-based system, 

2 comprising: 

3 computing a first plurality of wavelet and scaling coefficients from a plurality of 

4 first data points included in a first cartographic database; 

5 computing a second plurality of wavelet and scaling coefficients from a plurality 

6 of data points included in a second cartographic database; and 

7 generating the database error metric based on a wavelet transform involving the 
3 first and second pluralities of wavelet coefficients. 

3 

ffi 21 . The method of claim 20, wherein the error metric is a total error metric based on a 

's.Li! 

3 plurality of wavelet scales. 

W 22. The method of claim 20, further comprising: 

?4 selecting a wavelet scale; and 

fl restricting the error computation to the selected wavelet scale to generate a layer 

|3f error metric. 

5 

1 23. The method of claim 14, wherein the data points are selected from the group 

2 consisting of coordinate pairs and coordinate triples. 

3 

1 24. A system for generating a database error metric, comprising: 

2 a first cartographic database for storing a first plurality of data points; 

3 a second cartographic database for storing a second plurality of data points; and 

4 a processor, operatively coupled to the first and second cartographic databases, 

5 configured to compute a first plurality of wavelet and scaling coefficients and a second plurality 

6 of wavelet and scaling coefficients, respectively, from the first and second pluralities of data 

7 points, the processor generating a database error metric based on the first and second pluralities 

8 of wavelet and scaling coefficients. 
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1 25. The system of claim 23 , wherein the error metric is a total error metric based on a 

2 plurality of wavelet scales. 

3 

1 26. The system of claim 24, wherein the processor is configured to restrict the error 

2 computation to a selected wavelet scale to generate a layer error metric. 

3 

1 27. The system of claim 24, wherein the data points are selected from the group 

2 consisting of coordinate triples and coordinate pairs. 

3 
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METHOD AND SYSTEM FOR WAVELET-BASED REPRESENTATION 
AND USE OF CARTOGRAPHIC DATA 



5 ABSTRACT OF THE DISCLOSURE 

Cartographic data is represented using wavelet coefficients. Specifically, a cartographic 
shape point database consisting of shape points and nodes stored on a per geographic feature 
basis is converted to a cartographic wavelet database consisting of wavelet and scaling 
coefficients also stored on a per geographic feature basis. The shape of each geographic feature 
10 can be reconstructed from the stored wavelet and scaling coefficients for that geographic feature. 
■5*5 The wavelet transform organizes the information in a function according to scale, a concept 
^4 which is related to resolution. The wavelet-based representation facilitates operations such as 
o line simplification (generalization) and zooming. Multi-resolution error metrics are defined 

using the wavelet representation for assessing the accuracy of sampled cartographic databases 
Jtt with respect to a ground truth database. 
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